Thermal leptogenesis in a 5D split fermion scenario 

with bulk neutrinos 



Jukka Maalampi 

Department of Physics, University of Jyvdskyld, Finland 
and Helsinki Institute of Physics, Helsinki, Finland 

and 

Iiro Vilja and Heidi Virtanen 

Department of Physics and Astronomy, University of Turku, Finland 

August 3, 2010 

We study the thermal leptogenesis in a hybrid model, which combines the so called 
split fermion model and the bulk neutrino model defined in five dimensional spacetime. 
This model predicts the existence of a heavy neutrino pair nearly degenerate in mass, 
whose decays might generate a CP violation large enough for creating the baryon asym- 
metry of the universe through leptogenesis. We investigate numerically the constraints 
this sets on the parameters of the model such as the size of the compactified fifth dimen- 
sion. 

1 Introduction 

The origin of the baryon-antibaryon asymmetry observed in our Universe is one of the 
most intriguing open questions of the modern cosmology. It is also a question of particle 
physics as it is one of the most compelling pieces of evidence of the incompleteness of 
the Standard Model (SM). Had the early universe been matter-antimatter symmetric 
at the temperatures above the electroweak phase transition temperature (9(100 GeV), 
one would expect the ratio of the present number densities of matter over photons to be 
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nfl/n 7 ~ 10 18 (see eg [lj). This contradicts the observational result 

— = (6.1 ±0.3) x 1(T 10 (1.1) 

obtained by the Wilkinson Microwave Anisotropy Probe (WMAP) [2]. It has turned out 
that in order to explain this huge discrepancy one has to consider particle physics models 
that produce a larger asymmetry between matter and antimatter than what is possible 
to achieve within the SM. 

There are three general conditions [3], called the Sakharov conditions, that must be 
fulfilled for a baryon asymmetry to be created in the early universe: C and CP violation, 
baryon number (B) violation, and an exit from thermal equilibrium. Several particle 
physics models have been presented where the Sakharov conditions can be fulfilled with- 
out conflicts with other constraints, among these are various supersymmetric extensions 
of the SM ([3] and references therein). The SM is not among these viable models as 
the baryogenesis would require the Higgs boson to be much lighter than the experiments 
indicate [5]. 

One of the most appealing scenarios for the creation of the matter-antimatter asym- 
metry is the baryogenesis via leptogenesis (see eg [6J,[7J). In this scheme one extends 
the electroweak sector of the SM with interactions that violate lepton and baryon num- 
ber conservation. A net lepton number is generated perturbatively e.g. via triangle 
diagrams involving lepton number violating couplings, and this lepton number is then 
converted to a net baryon number by sphalerons through the so called Kuzmin-Rubakov- 
Shaposhnikov (KRS) mechanism [8J. Several models realizing this scenario have been 
proposed. Most of the contemporary leptogenesis scenarios, however, rely on a different 
mechanism, the so called thermal leptogenesis where a net lepton number is generated via 
heavy neutrino decays. In this mechanism, proposed in [6j , heavy neutrinos with a mass 
of the order of the Grand Unification scale (GUT scale) undergo CP violating decays 
producing an lepton-antilepton asymmetry among the decay products. Heavy neutrinos 
serve also another purpose in these models as they offer, via the see-saw mechanism [9J, 
an explanation for the lightness of the known SM neutrinos. 

Another class of models is provided by the so called low energy extra dimension brane 
models inspired by superstring theories. In these brane-world scenarios our universe 
is supposed to be a 4 dimensional hypersurface, called the brane, living in a larger 
dimensional space-time, called the bulk. It is supposed that the SM particles, including 
the ordinary left-handed neutrinos, reside on the brane, while sterile particles such as 
right-handed neutrinos are allowed to propagate also in the bulk |T0J [TT]. The extra 
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dimension theories can offer a new solution to the hierarchy problem by bringing the 
fundamental scale of gravity (M ) many orders of magnitude below the effective gravity 
scale, Planck scale Mpi. In addition, they might also explain the mass hierarchies within 
the SM fermion families [T2| IT3] . 

The flavour can be brought into extra dimension models for example by introducing 
a separate bulk neutrino for each SM neutrino [TU| [TTj . This scheme is not restrictive as 
far as neutrino mixing patterns are concerned allowing for a diversity of effective mixing 
matrices among ordinary neutrinos. Its shortcoming is, however, the great number of 
undetermined parameters it brings along, which makes its predicting power quite limited. 
Another possibility is to introduce just one type of bulk neutrino with flavour-universal 
couplings to the SM neutrinos [2], which means that the coupling between the brane 
and bulk neutrinos is the same irrespective of the SM flavour of the brane neutrino. It 
turns out, however, that this kind of scheme would lead to effective mixing matrices 
that are too rigid for reproducing the mass and mixing patterns of neutrinos observed 
in neutrino oscillation experiments. 

An extra dimension model that circumvents these problems was presented Dienes and 
Hossenfelder in [15J, where the bulk neutrino scheme is combined with the so called split- 
fermion scenario [13]. In the split-fermion scenario the SM fermions are each centered 
in the brane around a different locations and their mixings are due to the overlapping 
of their corresponding wave functions. The split-fermion scenario at such suffers serious 
fine tuning problems as the couplings between particles are exponentially sensitive to 
relative particle distances in the brane. In order the model to reproduce the observed 
features of neutrino mixing the relative locations of neutrinos on the brane are strictly 
constrained [16J. In the model proposed in [15] such fine tuning problems are avoided. 
The model is a hybrid model where the split-fermion picture is extended by including 
bulk neutrinos. It allows the effective neutrino mixing angles to be completely decoupled 
from the sizes of the wavefunction overlaps on the brane. 

In the present paper we will revive the hybrid model of [15] and study it from the 
point of view of leptogenesis. We will work with a simplified version of the model 
considering one extra dimension and just one neutrino flavour as flavour does not play 
any essential role in leptogenesis. Our aim is to investigate whether leptogenesis can 
be realized in the framework of the hybrid model, in particular whether the parameter 
values required by the leptogenesis scenario are in accordance with the general setup of 
the model. Finding the allowed size of the extra dimension is key and indicates whether 
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one dimension is sufficient from the beginning or not. This also reveals whether the 
model can address the hierarchy problem. The plan of the paper is as follows. In Section 
2 we will introduce the hybrid model in the form we shall use it. In Section we will 
consider neutrino phenomenology the model leads to. The realization of leptogenesis in 
the model is presented in Section 4. Section 5 gives a summary of the results and our 
conclusions. 

2 The Hybrid Model 

In this Section we describe the basic structure of the hybrid model following the original 
work [15| . The general framework consists of n/ neutrino flavour eigenstates = 
(u a , v a R) T (a — 1, . . . , n/) bound to live on the brane and a four-component fermion \I> = 
(ip + ,ij_) T that can propagate in the bulk. Apart from the four-dimensional spacetime 
there is one extra spatial dimension compactified with a radius R. 

It is assumed that each of the active brane neutrinos fy a has a coupling with the 
bulk neutrino \l/ through a Yukawa term a H Pr{^ + ^ c ), where H is a Higgs field. 
At this point we differ from the original model of (T5] by introducing a complex phase. 
We assume that there is a phase difference between the couplings of the the left-chiral 
ip+ and the right-chiral ip_ components of the bulk neutrino The phase is necessary 
for the leptogenesis as it allows for CP violation needed for the creation of a net lepton 
number in the decays of heavy neutrinos. 

Explicitly, the action on which we will base our analysis is written in terms of two- 
component spinors as follows: 

r n/ r 

J a=l 

+i>l{x,y)gh(x,y) [^j c + (x, y) + e iSa ^_(x, y)] } +h.c. 

Since we are dealing with a fat brane that is shifted away from orbifold fixed points, the 
coupling v\{x, y)(M*e lSa + gh(x, y)e lSa )ip-(x, y) is allowed [11] . Otherwise, if the brane 
was located at an orbifold fixed point, then the orbifold boundary conditions would 
forbid the brane neutrinos v a from coupling to the odd bulk modes of tp_(x,y). The 
universal coupling scale is M* = g(H) and the vacuum expectation value of the 5D 
Higgs is written as (H) = v/y/2irR, where g is the dimensionful Yukawa coupling in 
the five-dimensional spacetime and v is the vacuum expectation value of the Higgs field 
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in the ordinary four-dimensional spacetime. Hence the universal brane-bulk coupling 



strength has the expression M* = vgjyMR. 

For the gamma matrices T A = (F*, T 4 ) (jj, = 0, 3) we use the chiral representation 



H = 0, ...,3 



a' 4 ' 
W 

iIjm, 

The kinetic terms of neutrinos in the five- dimensional action are given by 



n j 

S„ = d 4 xdy^2vl(x,y)ia^d^ a (x,y), (2.2) 



S b = J d A xd y mT A d A y. 

The nonzero contribution to the kinetic term of the brane neutrino contains only the 
4-dimensional derivatives as the y derivative renders the integrand odd with respect 
to y and thus the y integral vanishes at its end points. We assume that the extra 
spatial dimension undergoes an orbifold compactification. By making use of the orbifold 
relations ip + (—y) = ip+(y) and ip__(—y) = —ip_(y) we can write the Kaluza-Klein (KK) 
expansions in the following form: 

*l>+(x,v) = -X=^f{x) + - 7 L= V^ n) (x)cos^, 



^Rfo " R 

For the brane neutrinos located in the fat brane we use the Gaussian wave functions 
u a (x,y) = -4^ ex P (~\ — 2°^ ) u »(x), (2.4) 



where v a {x) is a four- dimensional spinor. For simplicity we assume that the wave func- 
tions of all flavours in the brane have the same width of o <C R. We also follow the 
assumption that the Higgs field profile in the extra dimension is constant |13j and plug 
in the zero mode from the Kaluza-Klein expansion: 

h{x,y) = -j=h(x). (2.5) 

This choice ensures the canonical normalization of the kinetic term of the 4D Higgs field 
h(x). 
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3 Neutrinos in the hybrid model 

Let us study the neutrino sector of the hybrid model in more detail. Following the original 
analysis of [15] we determine the mass spectrum of neutrinos and the corresponding mass 
eigenstates. 

The mass matrix in four spacetime dimensions is obtained by integrating the actions 
S c , S u and Sb, given in Eqs. (12.11) and (12. 2p . over the extra dimension y leading to 

n f 



n>0 



a=l 

En 
R 

n>0 



hm 



n>0 



— ^\x) + y: ( 



V 



^X.t ,{n)c ( s , 

—ipl {X) + 



h{x)m% - {n) 



n>0 



(3.1) 
} +h.c, 



where nj is the number of flavours residing on the brane. For the volume-suppressed 
couplings between the fields on the brane and in the bulk we have used the following 
notations: 

gv 



m = M ± \ — — 



a 



a 



m 



mcos 



ny 



R 





n 2 cr 2 


j exp 


2irR 2 _ 



V2me zSa sin 



R 





n 2 cr 2 " 


j exp 


2irR 2 _ 



(3.2) 



In what follows we will assume that the brane-bulk coupling is weak and set mR <Cl. 
The mass terms appearing in the action (13. ip are collected together as to 



a=l 



S mass = / d A xl vl(x) m4 0)c + Y < + V>? )e (*) + m» _^\x) 



n>0 



n>0 



This can be presented in matrix form as follows: 



+ h.c. 
(3.3) 

(3.4) 
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where the mass matrix Ai is given by 

M ~- 



A^D(n/+l)> 



M 



L(n / +l)x(n / +l) 



/ ^ a 

I M L 


M D \ 
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•Mr/ 
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■ ■ ■ m L+ 




nt 
m l,+ 


n / 

77T.1 


■ ■ ■ "v+ 


rif 
m n[- 


V 
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m \ 








o r 







(3.5) 
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R(ooxoo) 
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o o I 

o o I 



V 
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The left- and right-handed fields are arranged into the vectors Al and Af£ as follows: 

K = 04 ^1° )C , ^+ )C , )C , -) T > 

■A/l = K,^ 0) ,^ 1) ,^ 1)C ,...,^ ) ,^ )C ,...) T 

The matrices .M and Al* can be transformed to a block-diagonal form by the trans- 
formation |15j 

-J 



T 



(3.6) 



where we have denoted 

The transformation takes the mass matrix Ai into the form 

'Ml 



M = T T MT 



M 



R , 



where 



Air 



Af L - kAId 

m a/ g m 
m 



En (K,- m n,+ + ?<,+ m n,-)- 



(3.7) 



(3.8) 



(3.9) 



rn 
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The sum over n in the left upper block m a p of the matrix M. L can be approximated 
by an integral over k = (a / R)n because the sum is rendered finite as the Gaussian width 
of the brane neutrinos a acts as a regulator. This results in 

Mia 



Mffiv« + vp))-M : £(v*-Vf>)) 



(3.10) 



+e 



Erf(^(^ + ^))+Erf(^(^-^)) 



where the error function emerges in the integrations over the trigonometric functions 
appearing in the quatities |15| . Upon block-diagonalizing Ai* the upper block 

becomes just the complex conjugate of .Ml- 

The transformation ( 13. 6 p renders the field vectors Ml and J\f£ to the form 



Ml 



{K^^ c ^\ ...,^ )c ,^ n) , -) T = T T K 



(3.H) 



where 



and 





= < + 


1 (0)c 




i ( n ) c 


/ ( n ) c 




= # 




= Vol- 


7(o) 


= Vi 0) , 



TP 



(«) 



(n)c 



(n)c 



R 



n 



a=l 

£ 

a=l 



a R c 

m n , + -v a 
n 



«-E«:-^ n) +<, + ^ n)c ) 



R 



n 



a=l 
n f 



R 



n 



a=l 



(3.12) 



(3.13) 



We proceed by determining the eigenvalues and eigenvectors of the neutrino mass 
matrix in the simplified case where we take into account the electron neutrino only 
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and assume that the muon and tau neutrinos are decoupled, in other words, the mixing 
angles between the electron neutrino and other active neutrinos in the brane are assumed 
to be zero. In order to make the scheme to work in the case of one extra dimension, 
two left-handed brane neutrinos are needed: in addition to the SM electron neutrino z/ e , 
there must exist a left-handed sterile neutrino. Introducing the additional brane neutrino 
ensures we get a spectrum where one mass eigenstate is light while the remaining two are 
heavy. In relation to sterile neutrinos, a recent neutrino experiment suggests there are 
no light sterile neutrinos [17] . In our case the extra brane neutrino is related to a heavy 
mass eigenstate. The requirement of a plausible mass spectrum in terms of leptogenesis 
compels us to set the number of neutrino flavours to two, that is rif = 2. 

With the above assumption, the mass matrix (13. 9p becomes a 3x3 matrix (rif = 2) 
of the form 



Mi 



( ran m i2 m ' 

m i2 m 2 2 m 
m m 



(3.14) 



y m m u J 

where m is defined in Eq. (13.21) and m a p in Eq. (l3.10p . The eigenvalues of .Ml are 

Ai ~ -(m u + m 22 ) - mi 2 , (3.15) 



A 2 
A3 



/-ll/ 
-V2m + -m 12 + -(mn + m 22/ 

r 1 1. 
V2m + 2 m i2 + ^( m n + m 22) 



and those of .M L are just A^, \\ and A3. Due to the complex nature of m a/ 3, these 
eigenvalues are generally complex. The physical masses are mi = |Ai|, m 2 = |A 2 | and 
m 3 — I A3 1 - One has mi « m 2 pa m^. The corresponding mass eigenstates are given by 
the following superpositions of the interaction eigenstates: 

1 i<? : 12 ( 1 m n ~ m 22 ~ c , 1 m n - m 22 ~ 

Xi = ~i=e 11 -uf + —=- -vl 

V2 V V2|mu-m 22 | V2|mn-m 22 | 

1 -a in ( 1 mTi — ma 9 _, 1 mti — m™ \ . 
V2 v V2|mn-m 2 2| V2|mn-m 22 | ' 

where 0j = arg(Aj). The complex factors in Xi,2,3 w iH give rise to the desired CP asym- 
metry in the decays of the heavy states \2 3- 
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4 CP Violation and Leptogenesis 



Let us move to study leptogenesis in the model described above. As was mentioned, an 
attractive scenario for creating the baryon asymmetry consists of generating a lepton 
number asymmetry through lepton number violating decays of heavy neutrinos, followed 
by the generation of baryon asymmetry from this lepton asymmetry via anomalous B+L 
conserving effects during the electroweak phase transition. 

The CP- violation in heavy neutrino decays arises in leading order through the interfer- 
ence of the tree level amplitude and the lowest order vertex corrections (6j [7J HHJ [191 120] - 
It has been shown, however, that in some cases the interference of the tree-level ampli- 
tude with the diagram where one heavy state is transformed to the other via light lepton 
and Higgs loop (called as the mixing amplitude), will give a major contribution to the 
CP violation |2Tj [22| |23| l2l] . This may happen if the states that mix are a pair of nearly 
degenerate heavy neutrinos. As we have seen, in the model we are interested in there is 
an almost degenarete neutrino pair X2, X3- We will therefore concentrate in what follows 
on the CP-violation arising from the interference of the tree level digram and the mixing 
diagram. 

Let us study the decays of \2- The Feynman diagrams relevant from the point of view 



of CP violation are those presented in Figures [TJ and HJ Fig 1 (a) depicts the tree level 



decay of \2 into a light neutrino and Higgs boson, and Figs 1 (b) and 1 (c) present one loop 



diagrams where the process proceeds through a transition of \2 into an intermediate %3- 
The corresponding diagrams for antineutrino production in the decays of \2 are presented 



in Figs. 2(a), 2(b) and 2(c) The diagrams for the decays of Xs are obtained from those 



of X2 decays, presented in the figures, by interchanging \2 and ^3. 

As the heavy neutrinos in the model we are looking at are nearly degenerate, their 
mass difference being |m 2 — m 3 | ~ M 2 a <C m, we can expect the interference between 



the tree level diagrams of Figs 1(a) and 2(a) and the one-loop mixing diagrams of Figs 



1(b) , 1(c) 2(b) and 2(c) to give the leading contribution to the CP violation in neutrino 



decays. 

The CP violation arises from the difference between the decay widths of the lepton 
and antilepton production channels. In terms of the amplitudes, the relevant quantity is 

\M + Mi\ 2 ~ \M + M~i\ 2 ~ 2Re(M*M 1 ) - 2Re(M* M 1 ) 

where Ai and Ai\ are the tree level amplitude and mixing amplitude, respectively 
and Aio and Aii denote the corresponding antiparticle amplitudes. The CP-asymmetry 
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parameter that takes into account the decay of both heavy mass eigenstates can be 
defined as 



r( X2 -> H^xil) - T(X2 -+ H X m) + T(x 3 gxnj - £( X3 -+ H X m) (4 ^ 

r(x 2 ->^txiL) + r(x2->^xiK) + r(x3->^ t xiL) + r(x3->^xiK)" 

(Another definition can be found in [18] where asymmetries are calculated separately for 
X2 and xs an d then added together. ) 



X2,P \ Xl (^1,^2), Pi 



_1_ 



(a) 



if 

X2,P , v Xa 1 

+ ' x 

Xl(^ C l,2) 

Xl(vi,2), Pi 



(c) 

Figure 1: The relevant Feynman diagrams for the process % 2 ~ > XxlH '. The tree level 




diagram due to the decay of X2 to a neutrino and Higgs is in Fig 1(a) Fig 1(b) and 1(c) 
depict the mixing diagrams due to the the decay of X2 to a light neutrino and Higgs. 
The mass insertion occurs prior and after the loop, respectively. 

The phenomenological constraints of the value of the CP-violation parameter e are 
obtained by relating it to the observed baryon asymmetry of the universe. The lepton 
number generated via heavy lepton decays is related to the CP violation parameter e 
through (see eg. [25J and [6J) 

Y L ~ k——. (4.2) 

9* 

Here the parameter k is a factor that describe the dilution of the lepton asymmetry due to 
various lepton number conserving and violating processes taking place in the primordial 
plasma. Its value is estimated to be k — 10~ 2 — 10 _1 [251 EH EZJ. The lepton number 
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/ 

H 

i H,p H 

X2,p ! \xs(i#) ! 
X ' « ' 

-X <w Xl(^l. 2 ) 

X2,P \ 

(a) (b) 

X2,P ,' \ X3 / 




Xl(^l,2) 

>Xl(^ C l,2),Pi 



(c) 

Figure 2: The relevant Feynman diagrams for the process \2 — >■ XirH- The tree level 



diagram due to the decay of %2 to a light antineutrino and Higgs is in Fig 2(a) Fig 2(b) 



and 2(c) depict the mixing diagrams due to the the decay of X2 to a light antineutrino 



and Higgs. The mass insertion occurs prior and after the loop, respectively. 

Yi created is partially transformed to baryon number Y# due to anomalous electroweak 
processes. The lepton number and the net baryon number are related through |28j 



— ~~~ — - Vi; = (4.3) 
s c s — 1 c s — 1 s 

where s is entropy, % and are the baryon and lepton number density, respectively. 
In our case the factor c s is 

c = 8 "/ + 4 = 2 1 (44) 
22n 7 + 13 57' v ' ' 

yielding Y# = —201^/37. Taking the observational value for the baryon asymmetry, as 
given in Eq. (11.11) . and the relation s = 7.04n 7 between the entropy s and the photon 
number density n 7 , we will arrive at the condition 

6.1 x io- 10 = !*£ = 7.04/s:^^-. (4.5) 

TLry C s 1 

Given the estimated values for the parameter k ~ 0.01 — 0.1, we the obtain the following 
order-of-magnitude estimation for the allowed values of the parameter e: 

-e ~ HT 7 - 10" 6 . (4.6) 
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Since we work in a regime where the SM particles have acquired masses via the elec- 
troweak symmetry breaking, the sphaleron transition is not as efficient as it would be 
in the symmetric phase. Thus f)4.6p corresponds to the highest amount of CP-violation 
possible to produce in the model. 

When calculating the amplitudes, we work in the on-shell renormalization scheme and 
so the real (dispersive) part of the mixing/ self-energy loop vanishes when the propagator 
mass coincides with the renormalized mass m 2)3 . The couplings between T>1 2 and z/ 12 , 
which are allowed by the structure of the effective theory (13. 9p , can be neglaected in the 
leading order as the Yukawa couplings between the brane neutrinos are very small, of 
the order vn a ^jv <^m/v. 



The mixing amplitude of Figs 1(b) and 1(c) is given by 

1 m m* 



X2 



'11 



m 



2a/2 v \m n - m 22 



■12 ^i(0 3 -0 2 )/2 



e i(03-0 x )/2\ x 



uiPru 2 



A* 32 ml + m 2 m 3 A 32 - iA 33 AX 2 ml 



m 



\A 3 3\ 2 ml 



2imiKeA 33 



where 



.4 



m 



:V2 



A 



2567T v 2 



m 



m n - m 22 
\mii - m 22 \ 
ma ~ m 22 
2567T v 2 Vlmn — m 22 \ 



)(' 



i(6»i/2+e 2 /2~6»3)j 



'( e i(*i-ftO + e 



(4.7) 



(4.8) 



The antilepton decay mixing amplitude (Figs 2(b) and 2(c)) is given by 

1 m mn — ^22 



■A 2 



2\[2 v \m n - m 22 



■(■ 



Xd2-6 3 )/2 _ J(9i-6 3 )/2 



X 



UlP L U 2 



Al 3 ml + m 2 m 3 A 23 



iA 33 A* 23 m 2 



m 2 2 -m 2 3 - \A 33 



2 m 2 — 2zm2ReA 33 ' 



where A 23 



(4.9) 



v4 33 . The corresponding mixing amplitudes due to the decay of \3 are 
1 m mX 



X:', 



'ii 



m 



22 



2a/2 v \m n - m 22 \ 



(1 



e i(e 2 -0i)/2j x 



uiPru 3 



^23 m l + m 2m 3 A 23 - imlA 22 A 23 

2imlReA 



1717, 



mi\A 22 \ 



(4.10) 



22 



1 m mn — "^22 
2y/2 v \m n - m 22 \ 



(l_ e ^l-fc)/2j x 



uiPlu 3 



m|A 32 + m 2 m 3 A 32 



imlA 22 A* 32 



ms 



\A 22 \ 2 mj 



2im 3 ReA 22 ' 



where A 22 = A. 



32- 
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We find the following lengthy expression for the CP-violation parameter e defined in 
(TO: 



~(m 2 + m 3 ) x <{ [{ml -m\ - \A 33 \ 2 m 2 2 ) 2 + 4m 2 t (Rev4 33 ) 2 ] m 2 x 



cos ■ 



3 — 8 2 



m \~ m \- |A 3 3| 2 m2)(m2m3(Imy4 3 3 - ImA 22 ) 

+m2(ReA 33 ReA22 + ImA 33 ImA 2 2 - ReA 33 ReA 33 - ImA 33 IniA 33 + ImA 22 - Im/1 33 )) 
+2m 2 ReA 33 (m 2 (-IniA 33 ReA 2 2 + ReA 33 Iim4 22 - ReA 22 + ReA 33 ) 
+m 2 m 3 (ReA 33 - ReA 22 )) 



+ sin • 



Oo — 3 



m 



A 33 | m 2 )(m 2 m 3 (ReA 22 + ReA 



33 / 



+m^(ImA 33 ReA 22 - ReA 33 ImA 22 + ReA 22 + ReA 33 )) 
+2m2ReA 33 (-m 2 m 3 (ImA 22 + Iim4 33 ) + m^ReA^Re^ 

+Imyl 33 ImA 22 + ReA 33 ReA 33 + ImA 33 Imy4 33 + ImA 22 + Iim4 33 )) 



(4.11) 



+ [{m\ 



m 



\A 22 \ 2 mlf + 4m^(ReA 22 ) 2 ] m 3 



x 



cos ■ 



6 2 — 8 3 



m 3 — m 2 — \ A 22 \ m 3 )(m 2 m 3 (Iniy4 22 — ImA 



33 / 



+m3(ReA 22 ReA 33 + ImA 22 IniA 33 - ReA 22 ReA 22 - ImA 22 IniA 22 + lmA 33 - ImA 22 )) 



+2m3ReA 22 (m 2 m 3 (ReA 22 - ReA 33 
+Rev4 22 ImA 33 - ReA 33 + ReA 22 )) 
9 3 — 9 2 



-ReA 33 lmA 9 



m 3K~ rttJ.fi33im.fi22 



+ sin • 



v m 3 - m 2 



A 22 | m 3 )(m 2 m 3 (ReA 22 + ReA 33 ) 



+m 3 (ImA 22 ReA 33 - ReA 2 2imA 33 + ReA 33 + ReA 22 )) 
+2m3ReA 2 2(-m 2 m 3 (ImA 33 + ImA 22 ) + m 2 (ReA22Re,4 33 

+ImA 2 2lni74 33 + ReA 2 2ReA 2 2 + Im^Im^ + IniA 33 + ImA 22 )) 



We have treated oth the mass difference \m 2 — 777.3 1 and the Yukawa coupling squared 

m 2 /v 2 as perturbation variables as they are roughly of the same order of magnitude and 

both small in comparison with the masses of the decaying neutrinos: 

,2 



m 



-m %3 ~ \m 2 - m 3 -C m 2 , m 3 . 



(4.12) 



The corrections of first order in these parameters are taken into account in (14. lip , higher 
order corrections are negligible. 
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A few comments concerning our result are in order. The expression which we have 
obtained for the CP violation parameter e is more complicated than, for example, the 
expressions obtained in minimal models [22] [23J [21] based on SO(10) GUT. This is 
mainly due to the fact that in the model we are considering, based on the assumption 
of an extra spatial dimension, the CP-violation is generated in an energy scale where all 
fermions have achieved their mass as a result of the electroweak symmetry breaking. As 
a consequence, both brane and bulk neutrinos contribute to the mass matrix. Hence, our 
light mass eigenstate involves both two brane neutrino states, which results in a more 
involved combination of the Yukawa couplings (Aij) in the mixing loop. For example, in 
|24| the factors fulfill = A*- { and thus terms O(Afj) cancel, whereas in our case A^s 
are not symmetric (A^ ^ A* t ) and the cancellation does not occur. 

In order to have a viable mechanism for the creation of the baryon asymmetry, all 
three Sakharov conditions have to be fulfilled. The third condition requires that the 
expansion rate of the universe, given by the Hubble parameter H(T), must be greater 
than the tree-level decay rate of any L-violating process. This condition will in our case 
set a constraint on the 5D Higgs vacuum expectation value v. 

The dominant L violating processes in the present model are the heavy neutrino decays 
considered above. The third Sakharov condition then requires that the decay rates obey 
the condition (heavy neutrino mass denoted by m^) 



-itree 



2H(T = m N ), (4.13) 



which guarantees that heavy neutrinos are out of equilibrium when they decay. The tree- 
level decay rate is easily calculated to be (see eg [25j for decay rates of heavy particles 
and CP asymmetry producation) 

r tree = ^fl-cos^iY (4.14) 
64tt v 2 V 2 J ' v ; 

and the Hubble rate at the decoupling of the heavy neutrinos is given by |25| 

2 

H(T = m N ) = 1.73V£t/-> (4.15) 

Mpi 

where g* is the effective number of degrees of freedom at the stage of heavy neutrino 
decoupling. One can take g„ ~ 100. The condition (I4.13P becomes 

m 2 mjy / Q\ — : x ' 



, < 3.46 x 6477^^(1-003^ — —J . (4.16) 
v 2 Mpi •> / 
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The CP-violating parameter e depends on parameters yi$ = ^Jtt / (2a)yi t 2, 5i,2, R and 
the scale of the light neutrino mass M^a which we take to be ~ 1 eV |29| . In Fig [3] we 
present the values of e as a function of the size R of the extra dimension for three sets of 
constant values of the parameters yi^ and 5\^- The solid curve in Fig [3] corresponds to 
the set y~i = 1.0, ij2 = 2.0, 8i = tt/12 and 82 = 2tt/3, the dashed curve to the set y~i = 1.0, 
y~2 — 2.0, Si = tt/12 and 82 = 47r/3, and the dotted curve to the set yi = 1.0, y~2 — 2.0, 
81 = tt/12 and 82 = tt/2. Comparison of the plot with the constraint (14. 6 j) implies that 



-edO" 6 ) 




R(TeV _1 ) 



2xl0~ 13 5xlO" 13 lxlO" 12 2xl0" 12 5xl0" 12 

Figure 3: The CP violation parameter e as a function of the size of the extra dimension R, 
with -ei(yi = 1.0,y 2 = 2.0A = tt/12,5 2 = 2tt/3), -e 2 {yi = 1.0,3/2 = 2.0,^ = tt/12,5 2 = 
4tt/3), -£3(^1 = 1.0,y 2 = 2.0,8i = n/12,5 2 = tt/2). 

the allowed size range for the extra dimension for the given set of parameter values is 
R ~ 2.0 x 10" 13 TeV _1 - 4.0 x 10~ 12 TeV~\ Varying the phase angles 8i, 2 to different 
quadrants shifts the allowed values of R, and the plausible range for R is roughly from 
lO-^TeV" 1 to 10- n TeV _1 . Hence, sufficient CP violation can be produced in this model 
when the size of the extra dimension is in the range from roughly the Planck scale to five 
orders of magnitude larger than the Planck length and the masses of the heavy neutrinos 
are of the order of ~ 1 TeV. 

In FigHJ the CP- violating parameter e is plotted clS clS cL function of the (normalized) 
brane location y 2 for three sets of values of the parameters 5 1-2 , yi and R. The solid curve 
corresponds to the set 81 = tt/12, 82 = 2tt/3, yi = 1.0, R = 10~ 13 TeV^ 1 , the dashed 
curve to the set 81 = tt/12, 82 = vr/3, yi = 1.0, R = 10~ n TeV -1 , and the dotted curve 
to the set 81 = tt/12, 82 = 7tt/6, yi = 1.0, R = 10~ n TeV -1 . We can see from this plot 
that the values of the brane location of the neutrino z/ 2 that lead to an acceptable vales 



17 



of e depend quite strongly on the values of the phase angles Si. 



-e(i(T 6 ) 
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Figure 4: The CP violation parameter e as a function of the brane location y 2 f° r Hi — 1-0, 
Si = tt/12, 5 2 = 2tt/3 and R = 10~ 13 TeV" 1 ) (solid curve); y x = 1.0, 5 X = tt/12, 5 2 = tt/3 
and i2 = KT 11 TeV" 1 (dashed curve); y x = 1.0, S t = tt/12, 5 2 = 7tt/6 and i? = 10" 11 
TeV -1 (dotted curve). 

Fig [5] presents the CP violation parameter e as a function of the phase angle Si for 
three sets of parameters. The solid curve corresponds to the set = 1.0, y 2 — 2.0, <5 2 = 
vr/2, R = W- n TeV" 1 ), the dashed curve to the set {y x = 1.0, y 2 = 2.0, S 2 = 4vr/3, R = 
10- 11 TeV" 1 ), the dotted curve to the set {y x = 1.0, y 2 = 2.0, S 2 = 3tt/4, R = 10" 13 
TeV^ 1 ). The allowed values of the brane locations yi j2 and phase angles Si, 2 are quite 
restricted overall as only highly limited intervals of yi j2 and S X)2 with any given R lead 
to e of correct order of magnitude and correct sign. However, letting 5 12 vary leads to a 
periodic pattern of the values of e that are acceptable. 

5 Conclusions 

We have investigated, from the point of view of the leptogenesis, a model with one extra 
spatial dimension. The model, originally presented in [T5j, combines the so called bulk 
neutrino model and the split neutrino model. In this hybrid model different neutrino 
flavours are assumed to be in separate locations in a thick four-dimensional brane and 
in bulk there reside sterile neutrinos that couple with these brane neutrinos. Our study 
shows that in this model a CP violation large enough for leptogenesis to work can be 
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Figure 5: The CP violation parameter e as a function of the phase angle 5\, when 
(yi = 1.0, y 2 = 2.0, 5 2 = vr/2 and R = lO" 11 TeV" 1 ) (solid), (y x = 1.0, y 2 = 2.0, 5 2 = 
4tt/3, i2 = 10- 11 TeV- 1 ) (dashed), {y x = 1.0, y 2 = 2.0, 5 2 = 3vr/4, R = 10~ 13 TeV" 1 ) 
(dotted). 

created through decays of heavy neutrinos. We found that the size of the extra dimen- 
sion should be in the ranee 10~ 16 TeV -1 to 10~ n TeV -1 in order to ensure the correct 
magnitude of CP violation. 

A few concluding remarks are in order. The leading contribution to the CP violation 
arises from the amplitudes where the tree level diagram interferes with a one-loop self- 
energy digram where there is a transition between two almost degenerate heavy neutrinos. 
The effect of such amplitudes on the CP violation has been earlier studied eg. in [24], 
where it was found that the tree level decay width removes the singularity that occurs 
when the two heavy neutrinos are degenerate. In our model the coupling structure is 
different from that of the model of [21], and in our case the cancellation of higher order 
Yukawa terms do not occur but, on the other hand, these terms are perturbatively small. 

Extending the analysis of our model to the cases of more than one extra dimensions 
could also be worthwhile because then the sizes of the extra dimensions could be larger 
than the one found here. Also, as noted in [15], adding more dimensions would essen- 
tially change the neutrino mass spectrum. Namely, the brane-bulk coupling m would be 
suppressed relative to the brane-brane couplings m a /3, which gives rif mass eigenstates 
with masses ~ M^cr and one exceedingly light sterile state with a mass of the order 
of magnitude of a/R 2 . This hierarchy is likely to lead to different constraints on the 
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parameters compared to the ones found in this paper. 

Finally, we have restricted the analysis to the electron neutrino only, but we expect 
the conclusions would be qualitatively similar in the case of other neutrino types. 

Acknowledgements We thank Tuomo Puumalainen for assistance at an early stage 
of this work. One of us (HV) expresses her gratitude to the Magnus Ehrnrooth Foun- 
dation and the Finnish Academy of Science and Letters (the Vaisala Fund) for financial 
support. 

References 

[1] M. S. Turner in Intersection Between Particle Physics and Cosmology, Jerusalem 
Winter School for Theoretical Physics, Vol. 1, p. 99, eds T. Piran and S. Weinberg 
(World Scientific, Singapore, 1986) 

[2] WMAP Collaboration, C. L. Bennett, et al, Astrophys. J. Suppl. 148:1 (2003) 
D. N. Spergel, et al, Astrophys. J. Suppl. 148:175 (2003) 

[3] A. D. Sakharov, Violation of CP Invariance, c Asymmetry, and Baryon Asymmetry 
of the Universe, JETP Lett 5:24-27 (1967) 

[4] Stephan J. Huber, Thomas Konstandin, Tomislav Prokopec, and Michael G. 
Schmidt, Electroweak Phase Transition and Baryogenesis in the nMSSM, Nucl. 
Phys. B 757:172-196 (2006). 

M. Carena, M. Quiros, A. Riotto, I. Vilja and C.E.M. Wagner, Electroweak baryo- 
genesis and low energy super 'symmetry , Nucl. Phys. B 503:387-404 (1997). 
P. Huet and A.E. Nelson, Electroweak baryogenesis in supersymmetric models, Phys. 
Rev. D53:45 78-4597 (1996). 

Junhai Kang, Paul Langacker, Tianjun Li, Tao Liu, Electroweak baryogenesis in a 
supersymmetric U(l)-prime model, Phys. Rev. Lett. 94:061801 (2005). 

[5] K. Kajantie, M. Laine, K. Rummukainen and M. E. Shaposhnikov, The Electroweak 
Phase Transition: A Non-Perturbative Analysis, Nucl. Phys. B 466, 189 (1996) 

[6] M. Fukugita and T. Yanagida, Baryogenesis Without Grand Unification Phys. Lett. 
B 174:45 (1986). 

[7] M. A. Luty, Baryogenesis via Leptogenesis, Phys. Rev. D45:455-465 (1992) 



20 



[8] V. A. Kuzmin, V. A. Rubakov and M. A. Shaposhnikov, On the Anomalous Elec- 
troweak Baryon Number Nonconservation in the Early Universe, Phys. Lett. B155, 
36 (1985) 

[9] P. Minkowski, Mu ->■ E Gamma At A Rate Of One Out Of 1- Billion Muon De- 
cays?,Phys. Lett. B 67, 421 (1977); M. Gell-Mann, P. Ramond and R. Slansky, 
Complex Spinors And Unified Theories, in Supergravity, edited by P. van Nieuwen- 
huizen and D. Freedman (North- Holland, Amsterdam, 1979), p. 315; T. Yanagida 
in Proceedings of the Workshop on Unified Theory and Baryon Number in the Uni- 
verse, edited by O. Sawada and A. Sugamoto (KEK, Tsukuba, Japan, 1979); S. 
L. Glashow, The future of elementary particle physics, in Proceedings of the 1979 
Carg'ese Summer Institute on Quarks and Leptons, edited by M. L'evy, J.-L. Bas- 
devant, D. Speiser, J.Weyers, R. Gastmans, and M. Jacob, Plenum Press, New 
York, 1980, pp. 687- 713; R. N. Mohapatra, G. Senjanovic, Neutrino Mass and 
Spontaneous Parity Nonconservation, Phys. Rev. Lett 44:912-915 (1980). 

[10] N. Arkani-Hamed, S. Dimopoulos, G.R. Dvali, J. March-Russell, Neutrino masses 
from large extra dimensions , Phys. Rev. D65:024032 (2001) 

[11] K. R. Dienes, E. Dudas, T. Gherghetta, Neutrino oscillations without neutrino 
masses or heavy mass scales: A Higher dimensional seesaw mechanism, Nucl. Phys. 
B557:25-59 (1999) 

[12] N. Arkani-Hamed, S. Dimopoulos, G. R. Dvali, The Hierarchy problem and new 
dimensions at a millimeter, Phys. Lett. B429:263-272 (1998) 

N. Arkani-Hamed, S. Dimopoulos, G. R. Dvali, Phenomenology, astrophysics and 
cosmology of theories with submillimeter dimensions and TeV scale quantum gravity, 
Phys. Rev. D59:086004 (1999) 

I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos, G. R. Dvali, New dimensions at a 
millimeter to a Fermi and superstrings at a TeV, Phys. Lett. B436:257-263 (1998) 

[13] N. Arkani-Hamed, M. Schmaltz, Hierarchies without symmetries from extra dimen- 
sions, Phys. Rev. D61, 033005 (2000) 

E. A. Mirabelli, M. Schmaltz, Yukawa hierarchies from split fermions in extra di- 
mensions, Phys. Rev. D61, 113011 (2000) 

[14] K. R. Dienes, I. Sarcevic, Neutrino flavor ascillations without flavor mixing angles 
Phys. Lett. B500:133-141 (2001) 



21 



[15] Keith R. Dienes, Sabine Hossenfelder, Hybrid model of neutrino masses and oscilla- 
tions: Bulk neutrinos in the split-fermion scenario, Phys. Rev. D74:065013 (2006) 

[16] G. Barenboim, G. C. Branco, A. de Gouvea and M. N. Rebelo, Neutrino masses 
and lepton flavor violation in thick brane scenarios Phys. Rev. D 64 (2001) 073005. 

[17] M. Maltoni, T. Schwetz, Sterile neutrino oscillations after first MiniBooNE results, 
Phys. Rev. D 76, 093005 (2007) 

[18] M. Flanz, E. A. Paschos, U. Sarkar Baryogenesis from a lepton asymmetric universe, 
Phys. Lett. B345:248-252 (1995) 

[19] L. Covi, E. Roulet, F. Vissani, CP violating decays in leptogenesis scenarios, Phys. 
Lett. B384:169-174 (1996) 

[20] W. Buchmiiller, M. Pliimacher, CP asymmetry in Majorana neutrino decays, Phys. 
Lett. B43L354-362 (1998) 

[21] J. Liu, G. Segre, Reexamination of generation of baryon and lepton asymmetries in 
the early Universe by heavy particle decay, Phys. Rev. D48:4609-4612 (1993) 

[22] M. Flanz, E. A. Paschos, U. sarkar, J. Weiss, Baryogenesis through mixing of heavy 
Majorana neutrinos, Phys. Lett. B389:693-699 (1996) 

[23] L. Covi, E. Roulet, Baryogenesis from mixed particle decays, Phys. Lett. B399:113- 
118 (1997) 

[24] A. Pilaftsis, CP violation and baryogenesis due to heavy Majorana neutrinos, Phys. 
Rev. D56:5431-5451 (1997) 

A. Pilaftsis, Heavy Majorana neutrinos and baryogenesis, Int. J. Mod. Phys. 
A14:1811-1858 (1999) 

[25] E.W. Kolb, M. S. Turner, The Early Universe, Addison- Wesley Publishing Company 
(1993) 

[26] W. Buchmiiller, T. Yanagida, Quark lepton mass hierarchies and the baryon asym- 
metry, Phys. Lett. B445:399-402 (1999) 



[27] M.-C. Chen, TASI 2006 Lectures on Leptogenesis, hep-ph/0703087 (2007) 



22 



[28] J. A. Harvey, M. S. Turner, Cosmological baryon and lepton number in the presence 
of electroweak fermion number violation, Phys. Rev. D42:3344-3349 (1990) 
R. N. Mohapatra, X. Zhang, QCD sphalerons at high temperature and baryogenesis 
at electroweak scale, Phys. Rev. D45:2699-2705 (1992) 

S. Yu Khlebnikov, M. E. Shaposhnikov, The Statistical Theory of Anomalous 
Fermion Number Nonconservation, Nucl. Phys. B308:885-912 (1988) 
M. Laine, M. E. Shaposhnikov, A remark on sphaleron erasure of baryon asymmetry, 
Phys. Rev. D61: 117302 (2000) 

[29] J. Bonn et al., The Mainz neutrino mass experiment, Nucl. Phys. B(Proc. Suppl.) 
91,273 (2001) 



